Tip effects in scanning tunnelling microscopy of atomic-scale magnetic structures. 
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The spin-polarized scanning tunnehing microscope (STM) can in principle resolve not only the 
electronic, but also the magnetic surface structure. We model recent STM measurements achieving 
magnetic resolution on the atomic scale by a first-principles method. It is shown that the signature of 
a specific magnetic or non-magnetic STM tip can unambiguously be identified. It is also established 
that the model of Tersoff and Hamann would yield an electronic as well as magnetic contrast of the 
surface which is well below STM resolution. 
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Currently, intensive research is undertaken to under- 
stand the complex behavior of spin-systems in thin fer- 
romagnetic (FM) and antiferromagnetic (AFM) metal 
films. The motivation behind this effort is the ever in- 
creasing demand for denser storage of information on 
magnetic disks, the main storage devices in the computer 
industry. While traditionally the emphasis was on long 
range spin-interactions, because these can be influenced 
by controlled growth of magnetic and non-magnetic mul- 
tilayers, short range interactions have only recently come 
into focus. In this range neighboring atoms with differ- 
ent magnetic properties can act as nanomagnets, their 
magnetic orientation driven by competing exchange in- 
teractions between the electrons. The field gained a con- 
siderable boost last year by the first demonstration of 
magnetism on the atomic scale [0. It is expected that 
the packing density of information on ultrathin magnetic 
films can ultimately be brought down to the very level 
of single atoms. We may therefore not exaggerate by 
saying that nanomagnetism is today the most promising 
field within the burgeoning area of nanotechnology. 

The spin-polarized (SP) scanning tunnelling micro- 
scope (STM), which was used to perform the first local 
measurements of antiferromagnetism on an atomic scale, 
can also serve as a tool to probe complex spin structures 
on surfaces, as recently demonstrated In this case 
a magnetic tip tunes in on the spin states of electrons 
on the sample surface and reveals the magnetic, rather 
than the electronic, structure of the sample. But the cur- 
rent theoretical model, employed to elucidate the mag- 
netic contrast, disagrees with experiments by one order 
of magnitude. In order to exploit the full potential of this 
method for quantitative measurements of magnetism on 
the atomic scale it is necessary to understand in detail 
how this contrast actually arises. 

We shall show in this Letter that the missing informa- 
tion is the chemical and electronic structure of the tip. To 
this end we simulate STM measurements from first prin- 
ciples, employing magnetic and and non-magnetic model 
tips and the STM code recently introduced (bSCAN [|). 
In particular we analyze, whether (i) the corrugation is 
really substantially enhanced if a tungsten terminated 



tip is used [|l|; and (ii) whether contamination of the tip 
with surface atoms might change the images. The first 
point seems important because extensive simulations of 
STM scans on an iron surface brought no evidence for an 
actual enhancement of the corrugation [^. In addition, 
it was shown on metal alloy surfaces that a tungsten ter- 
minated tip in fact diminishes rather than enhances the 
corrugation ||]. Furthermore, we could demonstrate that 
relaxation effects in the close distance range enhance the 
corrugation by about 10 to 15 pm ||]. Given this evidence 
the only route for a fuller understanding are detailed sim- 
ulations, which we present in this paper. 

The paper is structured as follows. First we present 
results of our calculation of the W(llO) surface with a 
single Mn over layer. The corrugation of the surface, if 
evaluated from the local density of states, will be shown 
to be very small and in fact below the limit of STM reso- 
lution. Then we demonstrate that due to the influence of 
the electronic properties of a tungsten tip the corrugation 
is reversed and enhanced by a factor of about five. This 
contrast reversal due to tungsten tips is quite common on 
magnetic surfaces It is the first time, however, that 
we also obtain a substantial contrast enhancement. Fi- 
nally, we present the simulated images of the W(110)Mn 
AFM surface for an Fe terminated and a Mn terminated 
tip. Here we show that the tip most likely used in one of 
the experiments was contaminated by Mn atoms. 

The electronic structure of the sample and the tip sur- 
face were calculated within density functional theory in a 
two step procedure, chosen for efficiency. First, we calcu- 
lated the groundstate interlayer distances of the surface 
and a W(100)X terminated tip, where X is the single 
apex atom. For this part of the calculation we used the 
Vienna ab initio simulation program (VASP) This 
groundstate geometry was then used in the full-potential 
augmented plane wave calculation of electronic surface 
properties (FLEUR, contrary to the version described in 
[|| our code does not include spin-orbit coupling). This 
full-potential method is still the most precise numerical 
tool for a complex magnetic system. For both steps we 
employed a spin-polarized generalized gradient approxi- 
mation to exchange and correlation energy according to 
Perdew et al. (PW91) @. 
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The STM currents and corrugation were computed 
with a perturbation approach based on Bardeen's inte- 
gral . The program computes the integral numerically 
by summing up all states within a given energy interval. 
The interval chosen in our simulations was 50meV, it is 
larger than the thermal broadening of states at the tem- 
perature of measurements, which Heinze et al give as 
16K 0]. The reason for this increase is that we sought to 
avoid an unrealistically selective evaluation, which might 
be distorted by numerical limitations. 

We started our simulations by computing the ground- 
state properties of the W(llO) surface with a FM or AFM 
Mn adlayer. This adlayer grows pseudomorphically on 
the tungsten surface [|l^. Initially, we determined the 
equilibrium position of the ion cores. The relaxations 
were performed with VASP. The energy cutoff of the op- 
timized pseudopotentials 0] for W and Mn was set to 
Ecut = 230eV. We used a Monkhorst-Pack grid of 
(8x8x1) /c-points for the Brillouin-zone integration. 
The system was simulated by a repeated slab of seven W 
layers and one Mn layers on both sides of the tungsten 
slab. The vacuum layer separating the repeated slabs 
was approximately nine A thick, this distance is in every 
case sufficient to guarantee that adjacent slabs are com- 
pletely decoupled. The lattice constant of the tungsten 
lattice had previously been computed by minimizing the 
crystal energy within the PW91 functional H. There, we 
obtained a slightly larger constant of 3.19 A compared 
to the experimental value of 3.16 A. The interlayer dis- 
tance between the tungsten surface and the Mn adlayer 
is smaller than the ideal bulk interlayer distance, we ob- 
tain an inward relaxation of 5.3 % (interlayer distance of 
2.13 A versus 2.25 A of the ideal lattice). We considered 
two configurations: the FM surface and a c(2 x 2) AFM 
surface. We could detect no surface buckling due to the 
magnetic properties of the surface atoms. In both cases, 
the FM as well as the AFM surface, we obtain the same 
vertical positions for the two Mn atoms. 

The equilibrium positions were the input for the cal- 
culation of electronic groundstate properties using a full 
potential code (FLEUR). In this calculation we used the 
same number of k-points (16), and the same GGA ex- 
change correlation potential (PW91). Heinze et al. H] 
obtained the AFM surface as the magnetic groundstate. 
In our calculation the energy difference between FM and 
AFM ordering of -5 meV/Mn is below the level of preci- 
sion even of full-potential methods. Therefore we would 
conclude that both states are equally probable. How- 
ever, our calculation did not include spin-orbit coupling. 
Spin-orbit coupling on this surface favors AFM ordering 
[|l| and thus shifts the balance slightly towards the AFM 
configuration. 

The induced moment in surface atoms is 0.24 /i^, it de- 
cays into the bulk and is close to zero for the central layer. 
The magnetic moment of the Mn atoms is 3.58 /i^, which 
is a consequence of the higher interatomic distance in our 
calculation due to gradient corrections for the exchange 
correlation functional. For the FM surface we obtain a 



magnetic moment of 3.62 /i^. To evaluate the density of 
states (DOS) in the Mn layer we used a triangular mesh 
of 36 k-points. The DOS of the two atoms in the AFM 
layer is shown in Fig. |l|. The majority and minority 
bands are below and above the Fermi level, respectively, 
the spin splitting in our calculation is about 4 eV (middle 
of majority to middle of minority band). Even though 
the exchange correlation and the lattice constant differ in 
our calculation, the DOS is virtually identical to the one 
obtained by Heinze et al. [Q. In case of spin-polarized 
surfaces the contributions from spin-up and spin-down 
states are computed separately. The spin of a tunnelling 
electron is conserved, which means that we do not con- 
sider the possibility of spin-fiips during tunnelling. 

Here we present only the results of our simulations for 
the AFM surface. Initially, we simulated a scan with a 
clean tungsten tip. The electronic structure of this tip 
is characterized by a high contribution of dxz and dyz 
states [|l5[. Furthermore, the DOS in the vacuum has its 
maximum at + 0.7 eV, and it is rather low at the Fermi 
level (0.42 states/eV). For these reasons a clean tip does 
not generally produce images of high contrast. In our 
simulation this is reflected by the corrugation height of 
2 pm (distance 4.5 A, - 3 mV sample bias and 0.5 nA 
tunnelling current), shown in Fig. |^ (a). The corrugation 
is lower than the experimental values by about 13 pm. 
Considering that the separation between tip and sample 
is already at the lower limit of stability [||, we attribute 
this difference to atomic relaxations of the system. The 
paramagnetic tip does not differentiate between the two 
Mn atoms; both are equal in the experiments and 
the simulations. A paramagnetic tip thus shows only 
the chemical structure of the surface. In this context 
it seems interesting to differentiate between effects of the 
surface electronic structure and effects due to the overlap 
between tip and sample wavefunctions. For this reason 
we have also plotted the local density of states (LDOS) 
at the same distance (see Fig. || (b)). The corrugation 
in this simulation is below the resolution of even a high- 
precision STM Neither the atomic positions, nor the 
magnetic structure of the surface as it appears in the 
scans can be accounted for by the electronic properties 
of the surface alone. 

The experiments with a ferromagnetic tip used a Fe- 
coated tungsten tip. We model the tip by a tungsten film 
with a Fe apex atom. Since the chemical nature of the 
apex atom is commonly the most important parameter 
in scans on metals [^, a single Fe atom seems sufficient. 
It is unclear, however, whether the tip was contaminated 
in the experiments by a contact with the Mn adlayer s. 
Therefore we also account for this possibility by simulat- 
ing scans with a Mn terminated tip. The results of our 
simulations are presented in Fig. ^. It should be noted 
that given plots refer to a median distance of 4.5 A. The 
current value was chosen accordingly and is smaller than 
for the tungsten tip (Fe-tip 0.1 nA, Mn-tip 0.2 nA). The 
reason for simulating scans at a specific height rather 
than a specific tunnelling current is the absolute value 
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of the current obtained in the measurements (34 nA). 
In a simulation this high current leads to a tip-sample 
separation well below the level, where the system is me- 
chanically stable j|] . For this reason we estimate the sig- 
natures of different tips at a given tip-sample separation, 
which we know to be at the lower end of tip stability. 

The simulations with magnetic tips in every case re- 
veal the magnetic structure of the surface. This is due 
to the changed coupling of spin-up and spin-down states 
of sample and tip. A detailed account of the electronic 
structure of Fe-terminated and Mn-terminated tips is 
given in [|l6[. Fig. || (a) shows the simulation with a 
Fe-terminated tip. While the qualitative features of the 
simulation are equal to the experimental results the 
corrugation height is much higher than the experimental 
value. Judging from results on other metal surfaces, this 
feature would be quite unusual. In fact, a corrugation 
this much higher than the simulated value has so far not 
been observed. For this reason we also simulated scans 
with a Mn-terminated tip. This simulation is presented 
in Fig H (b). The qualitative features remain the same, 
but the corrugation height (18 pm) is now in accordance 
with experimental values, if we assume a partial depo- 
larization of tunneling electrons due to the orientation of 
the magnetic axes. A similar assumption, with a similar 
result (polarization of tip states 40%, corrugation height 
7 pm), was made in the previous paper by Heinze et al. 
[y . We therefore conclude that the actual tip in the mea- 
surements was contaminated by Mn atoms of the surface. 

It might seem hard to understand, why the actual cur- 
rent in the experiments is so much higher than in the sim- 
ulations (34 nA versus 0.1 to 0.5 nA). In the simulations 
we obtained the current values by assuming that the dis- 
tance between tip and sample is at the lowest limit of me- 
chanical stability, or about 4.5 A(core-core distance) 
Generally, the tunnel current on a metal surface under 
these conditions amounts to 1-10 nA in simulations 
That it is only one tenth of this value in the present situ- 
ation is due to the low density of states at the Fermi level 
(see Fig. [l]). This low density of states yields a number 
density at the vacuum boundary of n(r) ^ 2 x 10~^A~^, 
which is about 10% of the number density encountered 
at the vacuum boundary of transition metals 0. Note 
that the same reduction due to the low density is en- 
countered by every method of STM simulation, whether 
it is based on scattering theory or perturbation theory, 
as long as it relies on density functional theory. 

Within groundstate density functional theory only two 
effects could lead to substantially changed current values: 
(i) The change of the magnetic splitting due to spin-orbit 
coupling, (ii) The change of the exponential decay due 
to image potentials. Both effects seem negligible in the 
present situation. The first, because the density of states 
is the same, whether spin orbit coupling is neglected, as 
in our calculation, or whether it is included, as in the 
previous calculation [|l[. The second does not change 
the number density at the vacuum boundary, which is 
already too low for any agreement with experimentally 



measured currents. So that we arrive at the conclusion 
that no current treatment of STM, based on first prin- 
ciples density functional theory, is suitable to reproduce 
the measured current values. 

It is conceivable that the difference is due to excita- 
tions or dynamic processes during tunneling. The energy 
range, for example, of tunneling electrons could not be 
limited to the Fermi level, if the system is excited. This 
question has been analyzed theoretically by Todorov, and 
it was found that the temperature gain due to a bias volt- 
age in the range of a few mV is about lOOK [|l^. For a 
substantial change the energy of the tunneling electrons 
would have to be about 2eV below the Fermi level. Since 
this energy shift corresponds to a temperature change of 
more than lOOOK, we can exclude this possibility. Fi- 
nally, it seems possible that phonon excitations due to 
interactions between tip and sample have an effect on 
the tunneling current. Here, we find that the change of 
the tunnel current due to phonon excitations is less than 
10% even for highly polar molecular bonds [|l9[. In sum- 
mary, we can exclude that any of these effects may alter 
the tunnel current by two orders of magnitude. There- 
fore the disagreement between experiments and simula- 
tions concerning the measured current values cannot be 
removed in any treatment relying on groundstate density 
functional theory. Since this theory has been so effective 
in the qualitative and quantitative modelling of experi- 
ments in condensed matter physics, we think that the dis- 
agreement provides a good reason for a careful revision of 
the experiments. Considering the acceptable agreement 
between measured and calculated corrugations and pend- 
ing further experimental results we think that the main 
contributions to the measured current should be inde- 
pendent of the lateral position of the STM tip, i.e. due 
to a current background not covered in our simulation. 
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FIG. 2. Corrugation of the surface from a constant current 
contour (a) and from the local density of states (LDOS, (b)). 
The corrugation vanishes for the LDOS. The corrugation with 
a paramagnetic tungsten tip is about 2 pm (tunneling condi- 
tions: - 3 mV bias and 0.5 nA tunnel current). Note that 
the highest protrusion of the current contour is at the hol- 
low site of the W(110)Mn surface. Ferromagnetic (full arrow) 
and antiferromagnetic (empty arrow) surface atoms are not 
discriminated. 
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FIG. 3. Corrugation of the surface from constant current 
contours and two different magnetic STM tips. The Fe ter- 
minated tip resolves the surface magnetic structure with a 
corrugation of 73 pm (a). The Mn terminated tip reveals an 
identical structure but a substantially lower corrugation of 18 
pm (b). In both cases the tip-sample separation is 4.5 A. 



FIG. 1. Density of states (DOS) integrated over the Mn 
surface atoms. The DOS of the ferromagnetic (top) and anti- 
ferromagnetic (bottom) Mn atom has the same structure, the 
polarization of spin states is reversed. 
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